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We test the models of vacuum energy interacting with cold dark matter and try to probe the 
possible deviation from the ACDM model using current observations. We focus on two specific 
models, Q = SfiHpA and Q = SfiHpc. The data combinations come from the Planck 2013 data, 
the baryon acoustic oscillations measurements, the type-la supernovae data, the Hubble constant 
measurement, the redshift space distortions data and the galaxy weak lensing data. For the Q = 
SflHpc model, we find that it can be tightly constrained by all the data combinations, while for 
the Q = SPHpA model, there still exist significant degeneracies between parameters. The tightest 
constraints for the coupling constant are /3 = —0.0261 q o |3 (for Q = ifjHpA) and fi = —0.00045 ± 
0.00069 (for Q = SfiHpc) at the Icr level. For all the fit results, we find that the null interaction 
/3 = 0 is always consistent with data. Our work completes the discussion on the interacting dark 
energy model in the recent Planck 2015 papers. Considering this work together with the Planck 
2015 results, it is believed that there is no evidence for the models beyond the standard ACDM 
model from the point of view of possible interaction. 


Modern cosmology encounters a great puzzle in ex¬ 
plaining the late-time cosmic acceleration, as a universe 
with barotropic or pressureless fluids, evolving accord¬ 
ing to the laws of general relativity, cannot ht the ob¬ 
servational data [1, 2]. This puzzle is phenomenologi¬ 
cally solved by introducing a new component with neg¬ 
ative pressure, generally called dark energy. However, 
the fundamental nature of dark energy is still unknown. 
At present, a cosmological constant A with equation of 
state w = — 1 is the simplest candidate of dark energy, 
which, however, is plagued with the so-called fine-tuning 
problem and coincidence problem [3, 4]. These theoreti¬ 
cal problems provide cosmologists with the motivation to 
consider some complex theories for explaining the cosmic 
acceleration, such as the dynamical dark energy models 
or the modihed gravity (MG) theories. In particular, in 
order to solve the coincidence problem, cosmologists have 
widely investigated the interacting dark energy (IDE) 
scenario [5-62] , where dark energy directly interacts with 
dark matter by exchanging energy and momentum. 

Regardless of the theoretical problems of A, the A cold 
dark matter (ACDM) model is still the most competitive 
cosmological model and is now even taken as a prototype 
of the standard model in cosmology, as it can excellently 
fit most of the current data with the least free parame¬ 
ters [63]. Nevertheless, other complex models mentioned 
above are also not excluded by the observations. Given 
this fact, we still need to continuously test the valid¬ 
ity of the ACDM model. Undoubtedly, any deviation 
from it, if confirmed by observations, would be a ma¬ 
jor breakthrough in cosmology. In fact, some degree of 
tension between different observations has been reported 
for the ACDM model since the Planck 2013 results were 
published [64]. This tension has attracted much interest 
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from cosmologists. Many works tried to probe or dis¬ 
prove these possible deviations from the standard ACDM 
model (see, e.g.. Refs. [65-69]). Particularly, in recent 
Planck 2015 results, the Planck Collaboration specifically 
discussed this topic in Ref. [70], where a series of dynam¬ 
ical dark energy and MG models were constrained by 
Planck 2015 data in combination with other astrophys- 
ical observations. The constraint results showed that a 
ACDM universe still resides in the allowed region of pa¬ 
rameter space at 2cr level for most of the data combina¬ 
tions. 

Although a number of models were studied in the 
Planck 2015 results, the IDE models were somehow ex¬ 
cluded (only a specific coupled quintessence model was 
discussed). In fact, compared with the dynamical dark 
energy and MG models, the IDE models can provide 
more features to fit the observations. Namely, the IDE 
models can not only affect the Universe’s background 
evolution but also directly change the growth history. 
In fact, they are able to both enhance and suppress the 
structure growth (see Ref. [47]), whereas, generically, the 
dynamical dark energy models only affect the background 
evolution (expansion history)^ and the MG models only 
tend to enhance the structure growth. We notice that a 
special IDE model with dark energy characterized by a 
scalar field (namely, a coupled quintessence model) was 
studied in the Planck 2015 results [70]. However, unlike 
a general IDE model where dark energy is phenomeno¬ 
logically described by some fluid, such a specific coupled 
quintessence model [m{(j)) = and V (/) = Vo(j)~°‘] 

always tends to enhance the structure growth; i.e., the ef- 


^ For a dynamical dark energy model with w a constant or varying 
slowly with time, it predicts a growth index, 7 = 0.55 -h 0.05[1 -h 
w(z = 1)], with small deviation from the standard ACDM value 
0.55 [71, 72]. Thus, dynamical dark energy models generally do 
not induce significant new feature to the history of growth of 
structure. 
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fective gravitational constant, GeS = G(1 + /3^) [70], is 
always larger than the standard gravitational constant 
G, no matter what sign of coupling constant (3 is cho¬ 
sen. Besides, due to the degeneracy between the coupling 
constant and the potential slope of the scalar field, the 
one-dimensional posterior distribution for the coupling 
constant somehow depends on the prior of the potential 
slope. (See Figs. 21 and 22 in Ref. [70].) From the above 
analysis, we believe that the discussion on the IDE model 
in Ref. [70] should be extended. 

Our aim is to test if there is any evidence for the de¬ 
viation from the standard (six-parameter) ACDM model 
from the point of view of possible interaction between 
dark sectors. To take careful aim at such a test, we 
only consider a one-parameter extension to ACDM in 
this paper. In other words, we focus on a special type of 
IDE model, i.e., the decaying vacuum energy model (vac¬ 
uum energy decays to cold dark matter, or vice versa) 
where dark energy is the vacuum energy with w = —1. 
In these models, the only extra parameter relative to 
ACDM comes from the coupling between vacuum energy 
and cold dark matter. Constraining this extra parame¬ 
ter could provide a test for the possible deviation from 
ACDM. The energy conservation equations for the den¬ 
sities of vacuum energy (pa) and cold dark matter (pc) 
satisfy 

PA = Q, (1) 

Pc = -3i7pc - Q, (2) 

where a dot denotes the derivative with respect to the 
cosmic time t, H is the Hubble parameter, and Q denotes 
the energy transfer rate. The decaying vacuum energy 
models have been widely studied in the literature; see, 
e.g.. Refs. [73-79]. 

Since we currently have no fundamental theory to de¬ 
termine the form of the energy transfer rate, we consider 
two phenomenological forms of Q, namely, Q = 3f3Hp\ 
and Q = SPHpc, where /3 denotes the coupling constant. 
It can be seen that such models only have one more free 
parameter (3 compared with the standard ACDM model, 
and /3 = 0 recovers the standard ACDM universe. Thus, 
by constraining the value of /3, we can clearly see whether 
any deviation from the ACDM model exists. Even with 
one extra free parameter, these IDE models can also 
change the background and perturbation evolutions of 
the dark sectors, simultaneously, which can be seen from 
the above equations and the following discussions. 

In the decaying vacuum energy models, if the cosmo¬ 
logical perturbations are considered, the vacuum energy 
should anyhow be disturbed by the perturbation of dark 
matter. Then the vacuum energy would no longer be a 
pure background but has perturbation, in principle. In 
return, the evolution of dark matter perturbation will 
also be affected by the perturbation of vacuum energy. 
The above scenario is embodied in the covariant con¬ 
servation laws, — Qj for J = A and c, where 

denotes the energy-momentum tensor and is 
the energy-momentum transfer vector. In this work, we 


choose = —Qc = Q'^c with the four-velocity of 
cold dark matter, so that there is no momentum transfer 
in the rest frame of cold dark matter. 

Now that the vacuum energy has perturbation, the 
crucial problem of how to handle this perturbation oc¬ 
curs. Note that we must take this problem seriously, 
since the calculation of the dark energy perturbation in 
a conventional way has induced several instabilities. A 
well-known example is that the perturbation of dark en¬ 
ergy will diverge when w crosses the phantom divide 
w = —1 [80-83]. In addition, in an IDE model the 
cosmological perturbations will blow up on the large 
scales at specific parameter values [52, 84, 85].^ These 
issues may result from the incorrect calculation of the 
pressure perturbation of dark energy, as pointed out in 
Refs. [86, 87]. Under such circumstances, we handle the 
perturbation of vacuum energy in this work based on the 
parametrized post-Friedmann (PPF) approach [88, 89]. 
The PPF framework applied to the IDE models, com¬ 
pleted in Ref. [86] for the first time, can successfully avoid 
the large-scale instability of IDE models and help us to 
probe the full parameter space of j3. For more infor¬ 
mation about the PPF approach and the calculation of 
the cosmological perturbations in the IDE models, see 
Refs. [86, 87]. 

To solve the background and perturbation equations 
for these two decaying vacuum energy models, we modi¬ 
fied the public PPF [88, 89] and Boltzmann code GAME [90]. 
We use the public Markov-chain Monte Carlo package 
CosmoMC [91] to explore the parameter space. The free 
parameter vector is {wf,, Wc, OmCi U ''^s, ln(10^°As), /3}, 
where ujb and uJc are the physical baryon and cold dark 
matter densities, respectively, rig and ln(10^°As) denote 
the spectral index and the amplitude of the primordial 
scalar perturbation power spectrum, respectively, 0 mc is 
the approximation to the angular size of the sound hori¬ 
zon at last-scattering time, and t is the optical depth to 
reionization. We set /? as a prior [—0.015, 0.015] for the 
Q = 3H(3pc model and [—0.15, 0.15] for the Q = 3PHp\ 
model; other free parameter priors are the same as those 
in Ref. [64]. In our calculations, we always assume two 
massless and one massive neutrino species with total 
mass = 0.06 eV. In addition, we also set f( = 0 

and cr = 0.4, in accordance with those in Refs. [86, 89], 
where and cr are two parameters used in the PPF ap¬ 
proach: is a function relating the momentum density 

of DE with that of the other components, and cr denotes 
a transition scale below which DE is smooth. (Note that 
and Cr can be calibrated in practice [87]. However, 
the choices of their values mentioned above are accurate 
enough for current observations [89].) 

For the observational data, we employ the Planck data 


^ The cosmological perturbations are stable only when < — 1 for 
Q oc pc model, and in > — 1 and /3 > 0 for Q oc pde model, if w 
is close to a constant. 
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in combination with other cosmological probes. Since the 
Planck 2015 likelihoods are not yet available, we use the 
data from the Planck 2013 release including the tem¬ 
perature power spectrum data from Planck [64] and the 
polarization power spectrum data from WMAP [92]. We 
shall use “Planck” to denote the above data combina¬ 
tion. Besides, we also use the Planck leasing data [93] as 
an additional option. For other observations, we use the 
following data sets: 

• BAO: the baryon acoustic oscillations measure¬ 
ments from SDSS at z = 0.15 [94], BOSS at 
0 = 0.32 and 0 = 0.57 [95], and 6dFGS at z = 0.106 
[96]. 

• SNIa: the type-la supernovae data from the Joint 
Light-curve Analysis (JLA) sample [97]. 

• Hq: the Hubble constant measurement Hq = 
(70.6±3.3)kms-iMpc-i [98]. 

• RSD: the redshift space distortions data from 
BOSS at z = 0.57 [99]. 

• WL: the galaxy weak leasing data from CFHTLenS 
[100]. All these data are utilized in the same way 
as Ref. [70]. Note that when we use the RSD 
data, we do not use the BAO data at z = 0.57. 
For simplicity, we also use “BSH” to denote the 
BAO-|-SNIa-|-i7o combination. 

Now we start to discuss our fit results. First, we con¬ 
strain the Q = ij3H model. The detailed fitting results 
are given in Table I. The one-dimensional marginalized 
posterior distribution curves for /3 are shown in the upper 
panel of Fig. 1, and the two-dimensional marginalized 
posterior distribution contours including /3 are shown 
in the upper panels of Fig. 2. For this model, the 
Planck data alone and the Planck-l-WL data combina¬ 
tion cannot provide good constraints on the coupling 
constant. The tightest constraint result comes from 
the Planck-|-WL-|-BAO/RSD data combination, which 
gives P = —0.026lo;o53 at la level. This result slightly 
favors a negative coupling constant corresponding to 
vacuum energy decaying into cold dark matter. The 
Planck-l-BAO/RSD data provide a similar constraint on 
/?. While the Planck-l-BSH data give /3 = 0.020i'ig;Q53 at 
the Icr level, showing that a positive coupling constant 
(corresponding to cold dark matter decaying into vacuum 
energy) is slightly favored. For all the data combinations, 
we find that there still exist significant degeneracies be¬ 
tween parameters in this model. 

Next, we constrain the Q = ^PHpc model. The de¬ 
tailed fitting results are given in Table II. The one¬ 
dimensional marginalized posterior distribution curves 
for P are shown in the lower panel of Fig. 1, and the 
two-dimensional marginalized posterior distribution con¬ 
tours including P are shown in the lower panels of Fig. 2. 
A remarkable feature of this model is that the coupling 
constant can be tightly constrained for all the data com¬ 
binations. Even using the Planck data alone, we can also 
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FIG. 1: The one-dimensional posterior distributions for p. 
The upper panel corresponds to the Q = SPHpA model and 
the lower panel the Q — SPHpc model. 


obtain a tight constraint result P = —0.0021 ± 0.0011 
at la level. This phenomenon can be easily understood. 
Since the interaction starts to work in the early matter- 
dominated epoch (due to Q ex pP), the time when cos¬ 
mic microwave background (CMB) is formed, it may 
induce significant effects on CMB. Another interesting 
phenomenon is that all the data combinations consis¬ 
tently favor a negative P, showing that vacuum energy 
decaying into cold dark matter is more supported. For 
this model, we find that the tightest constraint comes 
from the Planck-|-BSH data combination, which gives 
P = -0.00045 ± 0.00069 at Icr level. 

Finally, we test all the combinations with the addi¬ 
tional information from Planck lensing data. The con¬ 
straint results are slightly improved, as shown in Tables 
1 and 11. For a more complete analysis on the contribu¬ 
tion of the Planck lensing data to the constraints on the 
IDE models, see also, e.g., Ref. [62]. Looking at the con¬ 
straints with and without the Planck lensing data, we 
find that the null interaction (/3 = 0) is still consistent 
with current observations at the Icr level for most of the 
data combinations. Our results, together with the Planck 
2015 results, show that currently there is no evidence for 
the models beyond the standard ACDM model. 
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FIG. 2: The two-dimensional marginalized contours (68% and 95% CL) for the Q — 3(3HpA model (upper panel) and the 
Q — ZPHpc model (lower panel). 



TABLE I: The mean values and la errors for the main parameters of the Q = SfiHpA model. Note that A^min = Xmin, acdm ~ 
Xmini a-nd “-flensing” means that the Planck lensing data are combined. 
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TABLE II: The mean values and Icr errors for the main parameters of the Q — ZfiHpc model. Note that A^min = Xmin, acdm ~ 
Xmini and “-flensing” means that the Planck lensing data are combined. 
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We also wish to compare our analysis with previous 
ones. In our previous works, Refs. [86, 87], we have es¬ 
tablished the PPF framework for the IDE models and, as 
examples, we have made brief analyses for the interacting 
wCDM models (for simplicity, we call them the IwCDM 
models, where w is a constant) with Q = ZjSHp,. and 
Q = We found that, for the Q cc Pc model, 

using the geometric measurements can tightly constrain 
P [86]; whereas, for the Q oc pA model, using the geo¬ 
metric measurements alone cannot constrain /3 well, and 
the growth data from RSD can break the degeneracies 
and improve the constraints significantly [87]. For both 
IwCDM models, we found that the current observations 
prefer /3 < 0 at ~ 1.5a level. However, this conclusion is 
drawn in the models where two extra parameters (i.e., w 
and (3) relative to ACDM are introduced. Since the cur¬ 
rent observations can only constrain some combinations 
of /3 with other parameters—in particular, the degener¬ 
acy between /3 and w exists in current data—we actually 
cannot differentiate the effects from /3 and w. Concretely, 
according to the current observations, /3 is in positive cor¬ 
relation with w] i.e., to preserve the observational quan¬ 
tities, decreasing w can be compensated by decreasing 
(3 [86, 87]. Thus, the current data favor simultaneously 
re < — 1 and f3 < 0 [86, 87]. In order to avoid such a de¬ 
generacy and to test the simplest or minimal extension 
of ACDM (from the perspective of IDE), we have made 
an analysis for the interacting ACDM models (referred 
to as the lACDM models for simplicity, which are also 
called the decaying vacuum energy models in this paper) 
where cold dark matter is allowed to interact with the 
vacuum energy, without introducing any additional de¬ 
grees of freedom. Our analysis results show that, when 
w is hxed at —1, ,5 = 0 will be consistent with the cur¬ 
rent observations at the Icr level. We find that for the 
lACDM model with Q (x pA, /3 still shows strong de¬ 
generacies with other parameters, e.g., Dm, erg, and Hq. 
Thus, for this model, more accurate RSD measurements 
are needed in the future to help break the degeneracies 
(note that in this work, we only use one RSD point, i.e., 
the BOSS CMASS measurements of the three parameters 
Dy/rdrag, Fap, and /erg, evaluated at Zes = 0.57). 

Recently, in Ref. [79], the authors investigated the 
lACDM model with Q oc HpA, where the coupling 
strength /3 is set to be redshift dependent and is taken 
to be different values in four redshift bins in the analy¬ 
sis. They used the Planck 2013 data in combination with 
the SN Union2.1 sample and the RSD measurements (10 
points) to constrain the model and found that the null in¬ 
teraction (/3 = 0) is excluded at the 99% confidence level. 
One of the motivations of this paper is to check whether 
the null interaction is excluded for the Q (x HpA lACDM 
model. But our results show that the null interaction case 


of /3 = 0 is still consistent with the current observations 
at the Icr level (except for Planck-|-WL). The difference 
between our results and those of Ref. [79] may come 
from the following aspects: (i) We only used one RSD 
point (BOSS CMASS measurement at ^eff = 0.57), but 
Ref. [79] used 10 RSD points. We also included the 
Alcock-Paezynski (AP) effect [101] in the RSD measure¬ 
ment. Our utilization of the astrophysical observations is 
in exact accordance with the Planck Collaboration [70]. 
(ii) We handled the perturbation of the vacuum energy 
within the PPF framework [86, 87], but in Ref. [79] the 
perturbation of vacuum energy was neglected, (iii) We 
assumed that the coupling strength /3 is a constant, but 
Ref. [79] assumed /3{z) is a binned (stepwise-defined) 
function. It is of great interest to carefully test these 
aspects, and we leave such an analysis to a future work. 

Testing the validity of the standard ACDM model from 
observations is a vital task in current cosmology. Among 
the models beyond the ACDM model, the interacting 
dark energy model can provide more features to ht the 
observations. In this work, we consider the simplest or 
minimal extension of ACDM from the perspective of in¬ 
teracting dark energy, i.e., the decaying vacuum energy 
models (or the lACDM models). We have studied two 
decaying vacuum energy models with the energy trans¬ 
fer Q = 3(3HpA and Q = 2>j3Hpc, respectively. Both 
models only contain one more free parameter 13 com¬ 
pared with the ACDM model. By constraining f3, we 
can clearly see whether any deviation from the ACDM 
model exists. We used the PPF approach to calculate the 
perturbation of vacuum energy. Our data combinations 
come from the Planck data, the BAO measurements, the 
SNIa data, the Hq measurement, the RSD data and the 
WL data. For the Q = 3PHpc model, we found that it 
can be tightly constrained for all the data combinations, 
while for the Q = 3/3HpA model, there still exist sig¬ 
nificant degeneracies between parameters. The tightest 
constraints for /? are (3 = —0.026/q q 53 (for Q = 3j3HpA) 
and 13 = -0.00045 ± 0.00069 (for Q = 3l3Hpc) at the Icr 
level. For all the fit results, we did not find any evidence 
beyond the standard ACDM model. 
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